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ABSTRACT 

We analyse the phase-space structure of simulated thick discs that are the result of a 
significant merger between a disc galaxy and a satellite. Our main goal is to establish 
what would be the characteristic imprints of a merger origin for the Galactic thick 
disc. We find that the spatial distribution predicted for thick disc stars is asymmetric, 
seemingly in agreement with recent observations of the Milky Way thick disc. Near 
the Sun, the accreted stars are expected to rotate more slowly, to have broad veloc- 
ity distributions, and to occupy preferentially the wings of the line-of-sight velocity 
distributions. The majority of the stars in our model thick discs have low eccentricity 
orbits (in clear reference to the pre-existing heated disc) which gives rise to a char- 
acteristic (sinusoidal) pattern for their line of sight velocities as function of galactic 
longitude. The ^-component of the angular momentum of thick disc stars provides a 
clear discriminant between stars from the pre-existing disc and those from the satel- 
lite, particularly at large radii. These results are robust against the particular choices 
of initial conditions made in our simulations, and thus provide clean tests of the disc 
heating via a minor merger scenario for the formation of thick discs. 

Key words: methods: N-body simulations - galaxies: formation, kinematics and 
dynamics, structure - Galaxy: formation, kinematics and dynamics 



1 INTRODUCTION 



Thick dis cs appear to be a c o mmon feature of disc 
' galaxies (|Yoachim fc DalcantonI |2006| . and references 
] therein). They are believed to be snap- frozen relics of 
1 disc salaxy formation that took place at high-redshift 
eeman fc Bland-HawthornI [200^ 1. This is because their 
dominant stellar populations appear to be old, as indicated 
for example by the co lours of the envelopes of ed ge-on 
disc galaxies (see, e.g. iDalcanton fc BernsteinI |2002| . and 
follow-up papers). More direct evidence comes from local 
studies of the Milky Way's thic k disc, whose stars a re 
mostly older than 10-12 Gyr (e.g. lEdvardsson et al.lll993 ). 

In the case of the Milky Way, the thick disc has struc- 
tural, kinematic and chemical properties that in general are 
significan tly different compared to the other Galactic com- 
pone nts (|Bensbv et all l2005l : lJuric et all l2008l : IVeltz et all 
I2OO8I ). This implies that it is quite likely it ha s followed it s 
independent evolutionary path, as d iscussed in Wys3 (120041 ') 
and references therein (but see also llvezic et al.ll2008l ). 

As stated above, the Galactic thick disc provides us with 
a window into the high-redshift Universe. However unlike the 
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Galactic stellar halo, which is also ancient, the thick disc 
contains a non-negligible fraction of the total stellar mass 
of the Galaxy (between 6 and 15 per cent, e.g. iRobin et al.l 
I2OO3I : [ Juric et al] 120081 ). what enhances its importance as a 
tracer of the events that took place at early epochs. On the 
other hand, and like for the stellar halo, imprints of that 
early history may be present in the form of dynamical or 
chemical substructure. This information could, in principle, 
be retrieved relatively easily from large surveys of nearby 
stars implying that we can hope to directly test the various 
scenarios proposed for its formation. 

These various scenarios may be classified according to 
the relative importance of dissipative pro cesses. The collapse 
of a gas cloud with a large scale-height (|Eggen et al.l 19621: 
Burkert et all Il992l ). or intense star formati on (iKroupal 
2002 ) , perhaps triggered by gas rich mergers jBrook et al 



Z(){J^ ) , pernaps triggered by gas ricn mergers l|Drooit et al. 
2004 ; iBournaud et al.l |2007| ) are two dissipationally driven 



the vertical heatin g of a thin 
aLl liooi . 



models. On the other hand, 

disc during a merger event (e.g. lKazantzidis et al.ll2008l . and 
reference s therein); or the d irect accretion of satellites pro- 
posed bv lAbadi et al.l (|2003| ) are examples of (mostly) dissi- 
pationless processes. 

In the case of a merger origin of the Galactic thick 
disc, one may expect significant substructure to be present, 
particularly in the phase-space distribution of its stars. In 
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fact, evidence of such merger debris in our Galaxy has been 
mounting over the past decade. Examples are the substan- 
tial group of stars located a few kpc from the Galactic plane 
with kinematics intermediate betwee n the canonical thic k 
disc and the canonical stellar halo (|Gilmore et alj l2002h ; 
or the significant asymmetry in the distribution of thick 
disc stars i n the first Galactic quad r ant with respect t o 
the fourth (|Parker et al.l |2003| . |2004| : iLarsen et alTl2008h . 
In addition to substr ucture of the thin/thick disc su ch as 
the Arcturus stream (lEggenlll996l:|Navarro et al.ll2004 ). the 



Monoceros ring llYannv et al.ll2003f r the Canis Major dwarf 



l|Martin et al.l 20041 ). distinctive stellar groups in the solar 



vicinity with peculiar ages, metallicities and kinematics have 
also been discovered (|Helmi et al.ll200^ . All these observa- 
tions thus appear to support models in which the thick disc 
was fo rmed by acc retion and /or merger events. 

In lVillaloboTfc Holmi (2008, hereafter, Paper I) we pre- 
sented a series of simulations of minor mergers between a 
disc galaxy and a relatively massive satellite, with the aim 
of modelling the format ion of a thick disc (a scenario first 
explored thoroughly by iQuinn et al.l 1 19931 , and which has 
received much attention since then). The present work is in- 
tended as a follow-up study where we particularly focus on 
characterising the phase-space properties of these simulated 
thick discs. Our ultimate goal is to find sets of observables 
that would allow us to recover traces of the merger that 
may have lead to the formation of the Galactic thick disc. 
In principle, it should also be possible to distinguish stars 
from the intruder satellite from those of the heated disc in 
the final aftermath. 

We are also motivated to pursue such a study at this 
point in time by the surveys that are currently mapping the 
Milky Way galaxy and its components in great detail. In par- 
ticular, the spectroscopic RAVE and SEGUE/SDSS surveys 
are bot h providing large samples of stars with accurate kine- 
matics llZwitter et al.ll2008l : lAdelman-McCarthv et al.ll2008l : 
iLee et al.l |2008| ). Such surveys will enable a more precise 
characterisation of the Galactic thick disc, and should allow 
us to test the various formation scenarios discussed above. 
On a slightly longer timescale, the space astrometric mission 
Gaia will provide fu ll phase-space coordina tes for hundreds 
of millions of stars (jPerrvman et al.ll2001^ . from which we 
should be able to establish how the Galactic thick disc was 
assembled. 

The outline of this paper is the following. In Section 2 
we briefly describe the numerical simulations that are used 
in this study. In Section 3 we characterise the velocity distri- 
bution in local volumes resembling the Solar neighbourhood. 
In this section we put special emphasis on understanding the 
predicted distributions of heliocentric line-of-sight velocities 
since, in practice, these can be obtained with high accu- 
racy for large samples of stars, and are available even at the 
present time. Finally Section 4 summarises our conclusions. 



2 SIMULATIONS 
2.1 Description 

In iPaper ll we have performed a series of dissipationless N- 
body simulations of a single merger between a pre-existing 
disc galaxy and a satellite, in order to study the formation of 



thick discs in a c ontext of disc heating. Among the main re- 
sults presented in |Paper~3 , we find that such mergers are able 
to produce thick discs that are both structurally and kine- 
matically similar to those observed in the Milky Wa y and in 
external galaxies (see also lVelazguez &: Whitelll999l . and ref- 
erences therein) . Structurally, the simulated thick discs have 
larger scale-lengths compared to the initial disc and their 
scaleheights are 3-6 times larger with a clear dependence 
on the initial inclination of the decaying satellite. When 
compared to observations of the Galactic thick disc, the 
simulations seem to favour mergers with low/intermediate 
initial inclinations for the formation of this component, as 
sug gested by the rneasured az/uR ratio (as d etermined by, 
e.g., Soubiran et al .l2003l : IVallenarr et al ]|2006l ) and from the 
presen ce and amplitude of vertical gradients in the mean ro- 
tation (|Girard et al.ll2006l : llveiic et al.ll2008l ). 

These simulations are therefore a good starting point 
to study the detailed phase-space structure of the remnant 
system, and to establish if it is possible to dynamically dis- 
tinguish stars from the pre-existing disc from those of the 
accreted satellite. Furthermore, they may also be used to 
develop indicators to test the validity of this particular for- 
mation scenario for the Galactic thick disc. 

Below we summarise our suit e of num erical simulations, 
but refer the interested reader to iPaper ll for more details. 

We have explored the following initial configurations for 
the merger: (i) the structure and kinematics of the primary 
disc resemble those of a) the present Milky Way ("z=0" 
experiments); or b) a disc at redshi ft one ("z=l") accord- 
ing to the model of IMo et al.l (Il998l ). The latter represents 
a likely formation epoch of th e Galactic thick disc as sug- 
gested by the age of its stars l|Edvardsson et al] 1 19931 ): (ii) 
two stellar morphologies for the satellite: spherical or discy; 
(iii) two total (and stellar) mass ratios between the infalling 
satellite and the host galaxy (10% and 20%); and (iv) three 
initial orbital inclinations of the satellites, in both prograde 
and retrograde directions with respect to the rotation of the 
host disc. 

Both the host galaxy and the satellite are modelled self- 
consistently with both star and dark matter (DM) particles. 
The stellar component of the satellite has structure and kine- 
matics that are consistent with the observed fundamental 
plane of dE-f-dSphs galaxies. The satellite is launched far 
away from the host disc (35-50 times the host disc scale- 
length, depending on the experiment) and has orbital pa- 
rameters that are con sistent with c osmological studies of 
infalling substructure (|Bensonll2005l ). We will use here the 
simulations of 20% mass ratio between the satellite and the 
host galaxy. This results in a sample of 24 simulated thick 
discs. The stellar components of the satellite and of the host 
disc are modelled with 10^ particles, implying that the satel- 
lite's stars are over-represented in number (by a factor of 
five) in the merger remnant. 

As a reference. Table [T] summarises both the structural 
and the kinematical properties of the final thick discs of 
I Paper j . The kinematical properties listed here were mea- 
sured at R=2.4i?D for both "z=0" and "z=l" experiments, 
corresponding to ~11 kpc and ~5 kpc, respectively. Between 
brackets we also quote these values for the "z=l" case at 
R=3.6-Rd, i.e. ~8 kpc. In the rest of this paper we focus 
mainly on the "z—1" experiments, since the final stellar mass 
of the remnant system, ~ 1.4 x IO'^'^Mq, is comparable to 
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Table 1. Properties of final thick discs produced in lPapeT^ 
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85.7(64.8) 


76.4(53.9) 


55.4(47.8) 




56.7(54.8) 


53.9(44.8) 


48.2(28.3) 




36.8(20.9) 


41.7(32.5) 


41.8(42.6) 




96.1(117.3) 


107.8(120.7) 


122.7(133.7) 



- Scale-lengths (-Rd) and scale-heights (^o) in kpc. Radial (crp^), 
azimuthal (cr^) vertical (ctz) velocity dispersions and mean rota- 
tion (v^) in km s"-"^. 

- Kinematics are measured at 2.4Rd for "2=0" and at 2ARy) 
{3.6Ru) for "z=l" experiments. 

- Only experiments with prograde spherical satellites with masses 
20% of M^isc host E'l'e listed. 

that estimated for the Galactic thick disc (e.g. iRobin et al] 
l2007f ). 

For completeness, here we briefly remind the reader of 
the properties of the Milky Way thick disc. The measured 
scalelength of the Galactic thick disc is comparable to that 
of the thin disc, i.e. in the range 2.8- 4.5 kpc; and its expo- 
nential sc ale-height zo=700-2000 pc (|Larsen fc Humphreysl 
l2003l : lJuric et al. 2008). Kinematically, the velocity ellip- 
soid of the thick disc in the solar neigh bourhood is observed 
to be (gR,arf„ffz) ^('65,54,38) kin s~^ (jChiba fc Beersll200ll : 
IVallenari et al.l200g : rVehz et al.l200^. Thick disc stars have 
a rotational lag of ~30-50 km s~^ (|Chiba fc Beers! l2000l : 
IVeltz et al.ll2008l ). 

A close look at Table [1] shows that none of our sim- 
ulated thick discs reproduce exactly the properties of the 
thick disc of the Milky Way. An additional point to bear in 
mind is that the structure of the merger remnants is likely 
to change if a new thin disc is formed from freshly accreted 
gas, as generally envisioned in the models explored here (e.g. 
iFreeman fc Bland-Hawthornll2002l ). However, even if the de- 
tailed properties differ, we believe that studies of the dy- 
namical phase-space structure of our remnants should give 
us insight into what observables can be used for example, to 
distinguish in-situ stars from those that have been accreted. 



2.2 Definition of local volumes 

Our interest lies in determining the structure of phase-space 
of our merger remnants especially in small volumes which 
may resemble the Solar neighbourhood. This is motivated 
by our aim to eventually compare our predictions to obser- 
vations, and because typically we only have (access to) full 
phase-space information for relatively nearby stars. 

Given that the scalelengths of our simulated thick discs 
are ~ 35% smaller than what has been estimated for the 
Milky Way, it is not straightforward to decide which ra- 
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Figure 1. Coverage of the spherical volumes defined to anal- 
yse the phase-space structure of the simulated thick discs. The 
figures show the final spatial distribution (in a face-on view) of 
"2=1" experiments with spherical satellites on prograde orbits 
for three different initial orbital inclinations. Volumes are cen- 
tred at ilv=2.4i?D (~5 kpc) from the galactic centre and are 2 
kpc radius. Stars belonging to the heated disc (above) and to the 
satellites (below) are plotted separately to highlight the different 
spatial distributions. 



dius would correspond to the "Solar circle". This is why 
we explore the phase-space structure in volumes located at 
two different radii, namely 2.4i?D ('~5 kpc) and 3.6i?D (~8 
kpc) , where _Rd is the final thick-disc scalelength in our ex- 
periments (see Table [!}. 

To obtain phase-space information from particle sam- 
ples inside such volumes for our simulated thick disc, we 
first define a Cartesian system with the origin in the cen- 
tre of mass of the thick disc and the z-axis aligned with its 
rotational axis. The angular momentum vector points to- 
wards the galactic south pole {—z), hence the rotation is 
clockwise. We place four identical spherical volumes onto 
the disc plane defined as being perpendicular to the angular 
momentum vector. The centres of these spheres have coordi- 
nates (-fi?v,0,0), (-i?v,0,0), (0,+i?v,0) and (0,-i?v,0). -Rv is 
defined as either 2.4J?d or 3.6Rd- We explore spheres of sev- 
eral sizes but, unless noted otherwise, the rest of the paper 
will refer to volumes of 2 kpc radius. 

Fig.[T]illustrates the volumes' coverage of the final thick 
discs for the case of the "2=1" experiment using a prograde 
spherical satellite with initial inclination of 30° . The volumes 
are located at Rv=2.4Ro from the galactic centre and are 
labelled from 1 to 4 according to where they are placed, i.e., 
i+x,+y-x-y) = {l,2,3,4)- 

In the next Section we will characterise the velocity dis- 
tribution of stars in these local volumes. Because our rem- 
nant discs are not fully axisymmetric as evidenced in Fig. 
[1] it is important to establish first how the properties of the 
velocity distribution function depend on location. Therefore 
in Section [3T] we quantify the deviations from axisymmetry 
and establish their impact on the velocity distribution. We 
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Figure 2. Number of disc and satellite stars witliin each volume for all the experiments at "2=0" and "2=1". The volumes are labelled 
from 1 to 4 according to their location on the xy plane {+x,+y,—x,—y, see Fig. [TJ. In general, the spatial distribution of stars is not 
uniform for all 4 volumes, indicating different degrees of deviation from axisymmctry in both the heated disc and the satellite debris. 



then focus on the dynamical properties of the in-situ stars 
versus those accreted in Section [3.21 and Section [3.31 



3 RESULTS 

3.1 Effect of non-axisymmetries 

3.1.1 On the spatial distribution 

Fig. [1] shows that the distributions of heated disc particles 
and satellite particles are not symmetric with respect to the 
rotation axis of the system. Note that there does not seem 
to be a correlation in the deviations from axisymmetry be- 
tween the disc and the satellite stars. This suggests that 
these asymmetries are likely to have different origin. 

Fig. [5] gives a general overview, for all 24 of our exper- 
iments, of the distributions of stars in the final thick discs. 
This figure shows the variation of the number of stars from 
volume to volume for each experiment. Note that, in order 
to facilitate the comparison, the number of satellite stars has 
not been normaUsed according to mass ratio bet ween hos t 
disc and the stellar component of the satellite fsee lPaper j) . 
meaning that they are over-represented by a factor of five. 
Overall, in each experiment, both disc and satellite present 
some degree of volume-to-volume alternation in the num- 
ber of stars they contribute, demonstrating their asymmetric 
spatial distributions. In general, this alternation is different 
for the disc and for the satellite stars. 

Fig. m shows the global amplitudes of the m — 2 devi- 
ations of both disc (left panels) and satellite (right panels) 
spatial distributions as a function of the initial orbital in- 
clination of the satellites, for all 24 experiments. The am- 
plitudes are measured by first binning the final spatial dis- 
tributions of stars in cylindrical shells, within \z\ < 5 kpc 
and \z\ < 3 kpc for experiments at "z=0" and "z=l" , re- 
spectively. Then for each bin the Fourier components of the 
second harmonic of the angular distribution are computed 
as: 

1 ^ 1 ^ 

02 = ^ sin(2(?!)i), ^2 = ^ cos(2<^,), (1) 

i=l i = l 

l|Valenzuela fc Klvpin|[2003l and references therein) with 
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Figure 3. Global amplitudes of the m = 2 deviations from ax- 
isymmetry in the final spatial distributions of both disc (left pan- 
els) and satellite (right panels) stars for all experiments at "2:=0" 
(upper panels) and "2=1" (lower panels). 



the number of stars in each bin and the stars' angular po- 
sition. For each bin the amplitude of the second harmonic 
is A2 — (a2 -I- 62)72. The global amplitude is then defined 
as {A2) = -y/ {A2), averaging over radii R < 1.5Rd. Ad- 
ditionally, the angular phase of the m = 2 deviations are 
measured by finding the principal axes of the inertia ten- 
sors of the spatial distribution of stars projected onto the 
xy plane. In general, strong disc asymmetries are present 
out to a distance of ~3-4 kpc (~ 1.7Ru) for the "z=l" 
experiments. 

Figure |3] shows that in general, for both disc and satel- 
lite stars, there are no clear dependencies of the m = 2 am- 
plitude with inclination. However it is interesting to notice 
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Figure 4. Radial (left panel) and tangential (central panel) velocity dispersions in each volume as a function of the amplitude of the 
m = 2 deviations from axisymmetry in the spatial distribution of disc stars for all prograde "z=l" experiments. Velocity dispersions are 
shown considering both disc and satellite stars (solid circles) and disc stars only (open circles). The panel on the right shows the mean 
rotational velocity in each volume as a function of the angular separation between the volume and the major axis of the "bar" formed 
by the disc stars, for all "2:=1" experiments. 



that discy satellites on retrograde orbits have induced the 
largest m = 2 deviations on the heated discs. This can be ex- 
plained by the fact that this type of satellite loses mass faster 
as it evolve s in comparison to spherical satellites (see Fig. 3 
m iPaper j) . This results in a smaller drag force and hence 
in a longer decay timescale. Moreover, in general satellites 
on retrograd e orbits a re found to need more time to decay 
(see Fig. 2 in iPaperT) which is a consequence of the weaker 
dynamical friction exerted by the mai n disc as compared 
to th e prograde orbit case (see also IVelazguez fc White! 
1 19991 ). Therefore, the aforementioned satellite configuration 
induced a longer lasting and more powerful asymmetric per- 
turbation on the main disc, eventually leading to the forma- 
tion of a bar with a larger amplitude. 

On the other hand, the m = 2 deviations for the satellite 
debris seem to reach similar amplitude independently of the 
satellites' morphologies and initial orbital inclinations. This 
suggests that m = 2 deviations for disc and satellite stars 
may be caused by different dynamical effects. The satellites 
in our simulations have rather r adial orb its (apo-to-peri ra- 
tios ^ 5, see also Section 3.1 in IPaper t \. for all initial or- 
bital inclinations. This implies that their angular momenta 
are relatively low, or in other words that the amount of en- 
ergy associated to tangential motions is much smaller than 
that associated to motion in the radial direction. This in 
turn, means that the precession rate of these orbits is slow. 
Furthermore, our satellites are relatively cold in comparison 
to the host galaxy. These factors make the range of orbital 
angular frequencies of satellite stars quite small, implying 
that that the mixing in the angular direction has to take 
place on fairly long timescales. Thus as the satellite gets dis- 
rupted, its stars remain somewhat constrained in azimuth, 
which results in the observed m = 2 deviations, in rather 
close analo gy to what happens as a re sult of the radial orbit 
instability feinnev fc Tremaind [20081 ). 

There are two important implications of the charac- 
terisations made above. Firstly thick discs could well be 
asymmetric because of the presence of a small central bar. 
Secondly, the fraction of accreted stars may well vary with 
azimuthal angle in the plane of the thick disc. In recent 



years, iParker et al.] (|2003l . |2004 ) have found evidence that 
indeed the Galactic thick disc may be asymmetric. More re- 
cently, the suggestion has been made that the stellar asym- 
metry of faint thick-disc/inner-halo stars in the first quad- 
rant {I = 20° — 45° ) could be an in dication of a triaxia l thick 
disc or a merger remnant /stream (jLarsen et al.|[20o3) . 

3.1.2 On the kinematics 

In order to investigate the possible effects of the "bars" on 
the kinematics of stars in our local volumes, we have com- 
puted several statistical moments of the velocity distribu- 
tions (mean, dispersion, skewness and kurtosis). We have 
then correlated these properties to both the location and 
the amplitude of the m = 2 deviations of the spatial dis- 
tribution of disc stars (since this bar dominates in terms 
of mass/surface density as compared to the triaxial distri- 
bution of accreted stars which never do minate the central 
regions nor near the plane, see iPaper J) . In general, we do 
not find a significant dependence of the moments with re- 
spect to the strength of m = 2 deviations. This is shown 
explicitly in the left and central panels of Fig. |4j where we 
have plotted the radial and tangential velocity dispersions 
within each volume, for all prograde experiments at "z—V, 
as a function of the induced global m = 2 amplitudes. Note 
that similar results are obtained when the dispersions are 
computed considering either disc stars only (open circles) 
or disc and satellite stars jointly (solid circles). Only weak 
trends are found for the mean rotational velocity (v^) when 
computed as a function of the relative location of the volume 
with respect to the position angle of the bar, as shown in 
the rightmost panel of Fig. [3] When the various moments are 
computed in volumes which are located farther away from 
the centre, such as at 3.6i?D, the correlations are even less 
prominent. 

3.2 The kinematics of disc and satellite stars 

Fig. [5] shows the final velocity distributions of both disc 
(grey) and satellite (black) stars for the "2=1" experiments 
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Figure 5. Final velocity distribution of disc (grey) and satellite (black) stars for different initial orbital inclinations of a prograde 
spherical satellite: 0° (upper row), 30° (middle) and 60° (bottom); corresponding to "^=1" experiments. Stars are within spheres of 2 
kpc radius located at R=2.4i?j3 from the galactic centre. All satellite stellar particles have been plotted here, implying that they are 
over-represented by a factor of five. 




Figure 6. Final distribution of of disc (grey) and satellite 
stars (black; only one in five is shown) as a function of galactic 
radius, for prograde (top row) and retrograde (bottom row) satel- 
lites with initial orbital inclinations 0° (first column), 30° (second) 
and 60° (third). Stars are located within 2 < i? < 12 kpc and 
\z\ < 1 kpc. 



with a spherical satellite with initial inclinations of 0°, 30° 
and 60° on a prograde orbit. The stars plotted in this 



figure are enclosed within a spherical volume centred at 
R = 2ARu. 

This Figure shows that disc and satellite stars have radi- 
cally different distributions in velocity space. For instance, in 
the vb. — Vtf, plane, satellite stars are distributed in a banana- 
shape, whereas the disc stars define a centrally concentrated 
clump. This difference is a consequence of the more eccentric 
orbits of the accreted stars compared to those from the disc. 
The accreted stars that cross the volume under consideration 
have slightly different orbital phases, i.e. stars with wr = 
are at the apocentre while those travelling either towards (or 
away from) it have iir > (or hr < 0), leading to this ver y 
characteristic banana-shape (see also lHelmi fc Whitelll999l ). 
In the ur — plane satellite stars populate mostly the out- 
skirts of the distributions, with a clear dependence on the 
initial orbital inclination. Stars from the satellites on more 
inclined orbits show larger vertical velocities as they cross 
the disc plane as expected. Both the ur and distributions 
are very symmetric, showing that the stellar particles are 
well mixed in these directions by the end of the simulation. 
Satellite stars typically have lower mean rotational velocity 
in comparison to disc stars, as can be seen in the — 
plane. 

In Fig. |6] we show the final distributions of the z- 
component of the angular momentum (Lz) of both disc and 
satellite stars as a function of cylindrical radius, for "z=l" 
experiments with a discy satellite. In this figure, stars are 
located within 2 < i? < 12 kpc and |2| < 1 kpc. The trend 
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followed by the disc stars reflects a constant rotational veloc- 
ity (hence Lz cc R). The stars from the satellite, on the other 
hand, do not show this same trend, since their z-angular mo- 
menta are relatively constant with radius; the absolute value 
being related to the initial conditions of the satellite's orbit. 
This is why the separation between both types of stars has 
a mild trend with initial orbital inclination of the satellite, 
especially for prograde orbits, and also why the retrograde 
cases have negative (I/z). This figure suggests that in order 
to distinguish more easily between disc and satellite stars, 
it is better to study the kinematics at large galactocentric 
distances. 

Note that these conclusions can be generalised rather 
easily, and are quite independent of the particular initial 
configuration of the merger. For example, an increase in the 
initial tangential velocity of the satellite, will lead to a higher 
final value of Lz, implying that the distinction between ac- 
creted and in-situ stars will be more difficult near the Sun. 
In such cases, one would have to compile samples of stars 
located at much larger distances in the plane, such that the 
final Lz distribution at this location would be bimodal, as 
seen for example, at i? ~ 8 kpc in Figure [6] 

It is also important to stress here that the trend of Lz 
as function of R observed for the stars in the heated disc, 
constitutes a rather clean test of the formation scenario for 
the Galactic thick disc studied here. Such a trend would be 
expected only in the case a pre-existing thin disc was present 
at early times, since only then most of the stars would be 
moving on nearly circular orbits. 

Fig.[7]shows the Toomre diagrams for the "z=l" exper- 
iment with a prograde discy satellite with initial inclination 
30°, including disc and satellite stars. Note that the satel- 
lite particles are over-represented in this Figure by a factor 
five. The three velocity components have been convolved 
with Gaussian errors of 10, 20 and 30 km s"\ The volumes 
contain different number of particles due to the asymmetric 
spatial distributions discussed in Section 13.11 The separa- 
tion between both types of stars is evident even if the errors 
in the velocity are 20 km s~^ and should allow the distinc- 
tion between sa tellite stars and heate d disc stars in surveys 
such as RAVE l|Steinmetz et aLlbOOd ). 
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Figure 7. Toomre diagrams of "^=1" experiments for prograde 
discy satellites with initial orbital inclination 30° after convolu- 
tion of each velocity component with a Gaussian error of: 10 (up- 
per row), 20 (middle) and 30 (bottom) km s~^. In the case of er- 
rors '^10 km s~^, disc and satellite stars are easily distinguishable 
(disc stars have ~ 120 km s~^). Volumes 1 and 2 are shown 
to illustrate the difference in the distribution of stars due to the 
m = 2 deviations from axisymmetry in the spatial distributions 
of disc and satellite stars. Satellite particles are over-represented 
in this figure by a factor of five. 



3.3 Statistical tests on the velocity distributions 

3.3.1 KS test on the v-r, and distributions 

As just discussed, significant merger events like those simu- 
lated in |Paper~3 leave clear and long-lasting signatures in the 
kinematics of the final systems. Fig. [5] shows that disc and 
satellite stars clearly dominate different regions of velocity 
space, and that the final distributions are not Gaussian. In 
this section we quantify this non-Gaussianity, by studying 
the 1-dimensional velocity distributions along the R, (j> and 
z directions in our local volumes. To this end we perform 
the Kolmogorov- Smirnov/KuipeiEl statistical test (e.g., see 
iPress et al.lll992l ). Our null hypothesis is that the underly- 
ing velocity distributions are Gaussian along each one of the 
principal axes of the velocity ellipsoid {vr, v^, v^). 

^ The Kuiper variant of the test is chosen instead of the standard 
Kolmogorov-Smirnov since the first guarantees equal sensitivity 
at all values of the distribution, particularly in the wings. 



Fig. [S] shows the disc -f satellite velocity distributions of 
v-R and and the best Gaussian fit (dashed lines) for a 
"z—1" experiment with a prograde discy satellite with initial 
orbital inclination 30° . The stars considered here are located 
in volumes of 2 kpc radius at i? = 2.4i?D and R = 3.6i?D 
(left and right panels, respectively), and the number of satel- 
lite stars has been properly normalised (effectively, only 1 
in 5 stars is considered). This figure shows that in most 
cases our null hypothesis can be rejected with better than 
3(T confidence. This is mainly because the simulated velocity 
distributions have rather extended wings. These wings are 
preferentially dominated by satellite stars, also in the case 
of the distributions (see Fig. [SJ . However, in this latter 
case the contribution of satellite stars to the wings is too low 
(by number) to be detected by the test. As expected (Sec- 
tion the results of this test are nearly independent of 
the azimuthal location of the volumes. Note that even in 
these small samples, the features are very significant, and 
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Figure 8. Histograms of the radial (top) and tangential (bottom) velocities of stars within volumes centred at R = 2.4Hd and R = 3.6-Rdi 
for "2=1" experiments with a progradc discy satellite with initial orbital inclination 30°. KS/Kuiper tests have been performed in order 
to obtain a probability that the total distribution of stars (thick black line) is drawn from a Gaussian distribution (obtained as the best 
fit to the total distribution). Distributions of disc (grey line) and satellite (thin-black line) stars are also shown. The labels "KS Scr, 2a 
or ' X ' (meaning less than la)" indicate the confidence level with which one can establish that the distribution has not been drawn from 
the best fit Gaussian. Volumes 1 and 2 are shown to illustrate the difference in the distributions due to asymmetric spatial distributions 
of disc and satellite stars. 



should be relatively easy to detect observationally. Further- 
more, the velocity distributions do bear a clos e resemblance 
to tho se seen in the in-situ samples studied bv iGilmore et al] 
l|2002l ). 



3.3.2 Heliocentric line-of-sight velocities, vioa 

In practice, it is relatively difficult to obtain accurate full 
space velocities for large samples of stars. The main chal- 
lenge being to measure accurately their proper motions and 
distances for the stars. However, a possible workaround to 
this problem is to concentrate on heliocentric line-of-sight 
velocities («ios) which are comparatively cheaper to obtain 
and can be measured with high accuracy for large samples 
of stars. 

In our simulations we compute the vios as follows. First 
both the galactic longitude (/) and latitude (fe) are computed 
for each star with respect to the centre of the local volume, 
with the -|-a;-axis pointing towards the galactic centre and 
the -f j/-axis in the sense of rotation of the system. Then vios 
is computed as: 

"Uios ~ cos / cos b + Vy sin I cos b + sin b (2) 

where the three velocity components of each star are mea- 
sured with respect to the galactic centre (i.e. they include 
the galactic rotation). This implies that stars on circular 
low latitude orbits have V\ob ~ Vdrc sin / cos b, since Wx and 
should both be (close to) zero. 
The upper panels of Fig. [9] show the vios of disc (grey) 



and satellite (black) stars as a function of their galactic lon- 
gitude for the "z—V experiment with a prograde spherical 
satellite with initial orbital inclination of 30° . The most evi- 
dent feature is that stars of the heated disc are still rotating 
on fairly circular orbits in spite of the relatively massive mi- 
nor merger, as demonstrated by the sinusoidal dependence 
on galactic longitude of their line-of-sight velocities. 

We proceed to subtract the mean rotation of the sys- 
tem, Viot , in order to better quantify the contribution of the 
satellite stars to the wings of the distributions of wios. The 
mean rotation is found by fitting to the fios(0 distribution 
of both disc and satellite stars (after the proper normali- 
sation) the one-parameter function Ui.otsin(/). The result is 
the white curve in the top panel of Fig. |9l The middle and 
bottom panels of this Figure show both disc and stars radial 
velocities obtained after subtracting v^ot. The different be- 
haviour of disc and satellite stars is very clear now, and this 
is also the case for most of our experiments. In the particu- 
lar example plotted here, the maximum separation is found 
around I ~100-160° and I ~220-240° as illustrated in the 
figure. Note finally that the m = 2 deviations from axisym- 
metry do not significantly impact the distributions, i.e. there 
are no large volume-to-volume differences. 

In Fig. [To] we have plotted the vios distributions of 
disc+satellite stars, after subtracting Viot, for all "2=1" ex- 
periments with spherical satellites. The stars included in 
the histograms have galactic longitudes within an interval 
of SI — 40° around I = 140° (where the difference between 
the two populations of stars was largest according to the 
previous Figure). At these longitudes, I --^ 140°, most of the 
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Figure 9. Top panels: Heliocentric line-of-sight velocities of disc 
(grey) and satellite (black) stars. The white line shows a sin(i) fit 
to model the mean rotational velocity of the final thick disc. Mid- 
dle and bottom panels: The rotation signal has been subtracted 
for both disc and satellite stars, respectively. The figures corre- 
spond to the "2:=1" experiment using a prograde spherical satel- 
lite with initial orbital inclination 30° . Volumes 1 and 2 are shown 
to illustrate the difference in the distributions due to asymmetric 
spatial distributions of disc and satellite stars. In all the panels, 
the satellite stars are over-represented by a factor five. 



in the i-th bin. Therefore the numerator in Eq. ((3)l denotes 
the number of realisations for which the number of "random 
stars" is greater or equal the observed one for each bin. 

The resulting probabilities Pi{^ A^obs.i) are depicted as 
the solid circles in Fig. 1101 Typically very low probabilities 
(< 1%) are obtained in the wings of the velocity distribu- 
tion, as well as near the central peak. Recall that the wings 
are dominated by satellite stars (as indicated by the open 
circles). This implies that this test is both able to identify 
non- Gaussian features, as well as the presence of accreted 
stars. The low probabilities found near the central peak are 
due to the fact that this region is dominated by stars in 
a cold remnant disc, which survives the merger event, and 
contains 10 — 20% of the total mass of the thick disc (see 
Section 3.4 of lPaoer t i. 

Similar results are found for all our experiments (also 
those not shown in Fig. I10|l . In general, after analysing sim- 
ilar histograms for all the experiments, we find that even 
when the histograms contain a relatively small number of 
stars, it is possible to detect rare peaks in the wings which 
are composed mostly by satellite stars. 

At this point it is importa nt to rem ind the reader that 
the simulations performed in IPaper j do not include the 
growth of a fresh thin disc after the merger that led to the 
formation of the thick disc. This new, colder and (presum- 
ably) more massive thin disc would mostly dominate around 
wr ~ and Wz ~ km s~^. For example, in Figure |8] it 
should result in an enhancement of the central peak of the 
vk and distributions. On the other hand, in the direc- 
tion of rotation a pronounced peak should be present at 
i^ thin' where i^^t^in > '^.thick- This implies that the satel- 
lite stars are expected to remain the main contributors to 
the wings of all the velocity distributions. The deposition 
of a significant amount of mass in the galactic plane should 
lead to an increase in the rotational velocity of the thick 
disc stars, and hence the line-of-sight velocities should still 
show a sinusoidal dependence on longitude, but now with a 
larger amplitude than that visible in Fig. (9) Satellite stars 
will continue to dominate the wings of the vios distributions 
and should clearly become apparent after the mean rotation 
of disc has been subtracted. 



satellite stars have negative vios- This is manifested in most 
of the distributions shown in Fig. 1101 which are (slightly) 
asymmetric and tend to have a prominent negative velocity 
tail, particularly compared to the best fit Gaussian (denoted 
here by the dashed curves). To further emphasise this point, 
the open circles in Figure [TD] denote the fraction of satellite 
stars present at a given velocity (bin). 

We quantify the statistical significance of the features 
present in the vios distributions by computing a probability 
that measures the likelihood of the "observed" number of 
stars in a given velocity bin. We proceed as follows. We first 
generate A^tot.rcai ~ 10* random realisations based on the 
best Gaussian fit. We then define the probability of observ- 
ing Nohs,i stars or more in the i-th bin as: 

tot, real 

where Ng^cn^i denotes the number of "random stars" present 



4 SUMMARY AND CONCLUSIONS 

In this follow-up study we have analysed the phase-space 
pr operties of a sample of 24 simulated thick discs presented 
in iPaperll . These thick discs have been produced via a sig- 
nificant merger between a pre-existing disc galaxy and a 
satellite. In that study several combinations of the initial 
conditions of the progenitors were explored, such as: two 
redshifts of formation ("^=0" and "2=1"); two mass ratios 
between the main disc galaxy and the satellite (10% and 
20%); two different morphologies for the stellar component 
of the satellite (spherical and discy); and three initial or- 
bital inclination of the satellites (0°, 30° and 60°) in both 
prograde and retrograde orbits. 

The goal of this paper has been to find robust indica- 
tors of the merger origin of these simulated thick discs by 
characterising their phase-space structure. This involves es- 
tablishing which properties from the progenitors have been 
retained in the final system as well as ways to distinguish 
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Figure 10. Histograms of the heliocentric line-of-sight velocities after subtracting the mean rotation of the final thick disc, for "2:=1" 
experiments using prograde and retrograde spherical satellites with initial orbital inclinations of 0°, 30° and 60°. The histograms include 
disc and satellite stars within a slice around I ~140°, where the contribution of satellite stars to the wings is maximal (see Fig. |9]l. The 
width of Z-slices is 40° . The best Gaussian fits to the histograms are shown in dashed. The open circles show the fraction of satellite 
stars in each velocity bin. The solid circles denote the probability of the observed number of stars compared to what is expected from 
the best fit Gaussian at each velocity bin. 



dynamically in-situ and accreted stars. Our ultimate goal is 
to shed light onto the origin of the thick disc of the Milky 
Way by comparing the predictions of this model to already 
available (and future) surveys of nearby stars such as RAVE, 
SEGUE and Gaia. 

Our simulations show that the final spatial distributions 
of stars from both the heated disc and satellite are usually 
asymmetric with respect to the rotation axis of the system, 
but that each of these "m = 2" deviations are generally out 
of phase. The lack of apparent correlation between these 
asymmetries suggests that they could originate in difi'erent 
dynamical processes. Indeed, the "bar-like" distribution of 
heated disc stars seems to be induced by the asymmetric per- 
turbation of the decaying intruder, while the asymmetry in 
the distribution of satellite stars bears some analogy to the 
radial orbit instability. It is interesting to notice that simi- 
lar spatially asymmetric distributions have been observe d in 
the thick disc of the Milky Way (e.g. lLarsen et al.lbOOsI '). In 
terms of kinematics, these asymmetries are found to have a 
negligible effect on the velocity distributions of both heated 
disc and satellite stars. 

When samples of thick disc stars are selected in small 



volumes that resemble "solar neighbourhoods" , we find clear 
differences in the velocity distributions of in-situ and ac- 
creted stars. The stars from the heated disc are more cen- 
trally concentrated in ur and Uz while the accreted stars 
show broader velocity distributions in these directions. Ad- 
ditionally, the accreted stars rotate more slowly and they 
show a characteristic "banana-shaped" distribution on the 
■"R — Vif, plane. It is important to note that similar features 
arc observed independently of the size and location of the 
volumes on the midplane. 

The vertical component of the angular momentum, L^, 
as function of distance from the galaxy's centre is found to be 
a clear discriminator to separate heated disc stars from those 
accreted. In all our experiments in-situ stars have oc R 
whereas for accreted stars is approximately constant (i.e. 
independent of R) . This implies that the distribution of Lz is 
predicted to be bimodal, and that the bimodal nature should 
become more apparent with increasing distance from the 
galactic centre. The (x R behaviour thus provides a clean 
test of the presence of a pre-existing disc from which the 
Galactic thick disc stems. Note that such a behaviour would 
also be expected if the thick disc would have resulted from 
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re sonant interac t ions w ith transient spiral arms as proposed 
bv lRoskar et alj l|2008l ). In this case, however, one would not 
expect to find a second component associated to the accreted 
satellite. 

We also find that heliocentric line-of-sight velocities 
(viob) as a function of galactic longitude show that most 
heated disc stars remain on nearly circular orbits. This im- 
plies that even after a significant merger, the heated disc 
is able to retain part of the dynamical characteristics of the 
pre-existing disc (namely the relatively low orbital eccentric- 
ities). After subtracting the mean rotation, the wings of the 
Wios distributions are found to contain mostly accreted stars. 
Our analysis shows that the contribution of accreted stars 
to the wings is statistically measurable and, in principle, it 
could be easily detected by surveys in the solar neighbour- 
hood. 

Finally, it is important to highlight the robustness of the 
results presented in this paper regarding the different initial 
conditions explored. For instance, the significant separation 
between heated disc stars and satellite stars in terms of 
as a function of galactic radius. Even though the separation 
shows a mild trend with the satellite's initial orbital incli- 
nation, in all cases the separation is clear, especially when 
stars are located at large galactocentric radii. Similarly, the 
sinusoidal behaviour of vios{l) is observed in all experiments 
independently of the size or location of the volumes. The ro- 
bustness of these results allow us to consider them as direct 
probes of the disc heating scenario for the formation of thick 
discs. We aim to soon test these predictions on surveys of 
nearby stars such as RAVE to shed light on the formation 
of the Galactic thick disc. 
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